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Resume. Nous avons etudie la propagation d'un front de flamme dans une configuration unidimensionnelle en 
solvant numeriquement une equation de advection-diffusion- reaction (ADR) . Nous avons choisi un modele simplifie 
dans lequel le couplage avec d'autres phenomenes est neglige et le fiuide reactant est une mixture de gaz. En 
outre, la vitesse de I'ecoulement compressible est donnee de facjon analytique. A I'aide de la simulation numerique 
nous analysons la difference entre des modeles qui sont communement utilises dans les etudes fondamentales en 
combustion et en biologic. Puis nous avons caracterise I'effet de la compressibilite sur la dynamique du front et 
dans quelles conditions il est possible que la flamme s'eteigne avant la saturation de la reaction. 

Abstract. We have studied the front propagation in a one dimensional case of combustion by solving numerically 
an advection-reaction-diffusion equation. The physical model is siplified so that no coupling phenomena are 
considered and the reacting fluid is a binary mixture of gases. The compressible flow field is given analitycally. 
We analyse the differences between popular models used in fundamental studies of compressible combustion and 
biological problems. Then, we investigate the effects of compressibility on the front interface dynamics for different 
reaction types and we characterise the conditioins for which the reaction stops before its completion. 

1 Introduction 

Transport of reacting species advected by laminar or turbulent flows described by an advection reaction 
diffusion equation (ARD), is an issue of interest in many fields, e.g., population dynamics, propagation of 
plankton in oceanic currents, reacting chemicals in the atmosphere, ozone dynamics, complex chemical 
reactions, and combustion. While these phenomena have been and are being widely studied in the case 
of incompressible flows, the transport of reactive species advected by compressible flows is less discussed. 
Only recently has been extended to the study of population dynamics [1]. Combustion processes are 
particularly complex because they involve a large number of chemical species and a large number of 
reactions. In addition, the concentration of these species, and the release of chemical energy during 
reactions, greatly affect the flow field of the fuel and oxidizer mixture. Consequently, the numerical 
study of this phenomenon requires the solution of a large number of coupled partial difference equations 
(conservation of species, momentum conservation and energy conservation). 

The complexity can be greatly reduced if we minimize the physical details. In particular, if we neglect the 
effects of the coupling between the various equations and we minimize the number of species, combustion 
can be described approximately by an advection reaction diffusion equation: 

^ + V • (pu0) - V • (7V(/.) + s (1) 

which describes the spatio-temporal behavior of a fraction of the reactive mixture which moves with 
velocity u. The mean density of the mixture is p, 7 = pD is the diffusion coefficient and s is the source 
term. 

Although there is a strong connection between the study of population dynamics and combustion, as 
the simplified mathematical model is equivalent when we consider incompressible flows, in the case of 
compressible flows we will show briefly that the parallelism requires a little more attention. 
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2 Model 

If we consider a compressible reactive flow, and if the fluid is a mixture of N perfect gases, the local total 
mass balance over a volume V gives: 

|^ + V-(pu) = (2) 

where p{x, t) = pk is the density of the mixture and u(x, t) is the resultant velocity field. [2] [3] This 
equation is the summation over the N mass species conservation: 

^+V-[p(u + Vfc)cfc] = i?fe , k^l,...,N (3) 

where Ck{x,t) is the mass fraction (pk/p), Vfe(x, t) and Rk arc the diffusion velocity and the reaction rate 
of the A:*'* species respectively. Of course, by definition, J2k=i = J2k=i '^kCk = as mass cannot be 
generated during chemical reactions. 

If the mixture contains only two species, the pressure gradients are small, and volume forces are neglected 
[2] equation (3) for the fc*'' species can be expressed as: 



dck 



dck , ^ 
— + u • Vcfe 

ot 



V-{pDVck)+Rk , k = l,2. (4) 



where the diffusion velocity is rewritten according to Pick's law. 

Equations (4) describe the spatial, and time, behavior of the fractions (ci,C2) of a binary reacting mixture 
in a compressible velocity field. Such as, for example, the fresh air as well as the burnt gases in a 
combustion. 

If D is constant, equation (4) in conservative form becomes: 

■{puck)=D\/ ■{pVck)+Rk , k = l,2. (5) 

If D and p are constant, equation (4) become the widely studied incompressible advection- reaction- 
diffusion equation with constant diffusion coefficient. 
Equations (4) and (5) differ from the equation: 

^ + ^.{ue) = DoW^e + F{e) (6) 

which is often find in literature [1][4] to describe the advection, diffusion and reaction of a scalar in a 
compressible flow. This equation is typical in the study of population dynamics and the scalar 9{x. f.), 
is the concentration of a population [1]. However, this model is not correct for the concentration of the 
combustion products. In fact, in equation (6), if V-u 7^ the concentration can take values greater then 
one because is not a fractional parameter. Notably if one considers the Fishcr-Kolmogorov-Pctrovskii- 
Piskunov reaction rate (PKPP)[5], the rate of grow {F{9) = ac{l — c)) in equation (6), if not correctly 
rescaled, can take negative values. This could be in agree with the population dynamics point of view 
but is formally incorrect for the chemical reactions during a combustion. 



2.1 Problem statement and numerical details 

We have considered for simplicity a one dimensional case in which there is no feedback from the reaction 
rate, to the velocity field of the mixture. Moreover, the mixture is sufficiently dilute to neglect the heat 
release and the variation of temperature, thus we do not need to consider separately the energy equation. 
The velocity field is steady state and analytical so that the total mass balance becomes: 

= constant = /3 = 1 (7) 
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Figurel. Density (a) and velocity field (b) of the mixture as function of the parameter Uo- In figure (c) the three 
rate of reaction considered, as function of the mass fraction: the FKPP reaction (/i); the ignition type reaction 
(/2) with threshold Cs = 0.6 and the general mth order Fisher reaction (/s) with m = 10. 



where Uq is a parameter defined on the interval [0. 1] and L = 1 is the length scale of the velocity field. 
As shown in Fig.l, when C/q = the mixture moves with constant velocity, the density is constant and 
the problem becomes incompressible. On the other hand if Uq tends to one the density varies very rapidly 
between 0.5 and infinity. 

If wc consider the product pD = 7 constant and we collect and replace the velocity u with its analytic 
function, equation (4) becomes: 



dc 
di 



7 d^c 



dc 
dx 



(8) 



In this case, and later on, we will refer to c{x,t) as the mass fraction of the reacted products, like the 

mass fraction of burnt gases in a combustion reaction. 

If only D is constant the balance equation for the burnt gases (5) gives: 



dc+ 

~dt 



= u 



D d 




dc+' 


Jdx 




dx 



If we use equation (6) instead, we will get the following one dimensional ARD equation: 

^-D—-u—-9— fie) 
dt dx'^ dx dx 



(9) 



(10) 



We will consider three different non- linear rate of reaction, as shown by Fig.l.c 
The FKPP non linear rate of reaction or autocatalytic reaction: 

/i(c)=aic(l-c) (11) 

The ignition type rate of reaction: 

^^^""^ ^ { a2{c - -c)'/f 1>Z ^^^^ 
And the general m*'' order Fisher's non-linearity: 

/3(c) = a3c"(l-c) (13) 

If we use equations (8) and (9) c is defined in the interval [0, 1]. The reaction rates then, by definition, 
are always positives or zero which means that the reaction is irreversible. As we will show soon this is 
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not true for equation (6). 

The constants ai, 0.2 and arc chosen so that max {fi{c)) = ai/4. In this way reactions have different 
reaction rates but with comparable characteristic times. 

Equation (8), (9) and (10) have been solved numerically using a fourth order finite different scheme 

for spatial discretization and a fourth order Runge-Kutta method for time advancement. The numerical 
domain has a dimension = nL, where n is an even natural number, and periodic boundary condition 
have been applied. 



3 Results 

First of all we are interested in how the compressibility affects the time and the spatial behaviour of the 
fraction of burnt gas. We will solve now equation (8) setting up an initial condition for c{x, t) as: 

(14) 

Results are shown in Fig. 2. If we focus just the front interface, it is possible to identify a general behavior: 
when Uq = 0, the test case collapses to an incompressible advection-diffusion-reaction case where u{x, t) = 
1 and Do = ^ = 7 (as /3 = 1). Therefore, the burnt fraction moves with a constant velocity while it 
diffuses. The front interface is smooth end follows an exponential law. 
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Figure2. In contiunous line, the mass fraction of burnt gases after a transient of t = 15 in the case of FKPP (a), 
ignition (b) and generalized mth order Fisher reaction (c). In dashed line the initial condition where xo = 20.5 
and (To = 10. The parameter ^ = 0.1 (/3 = 1,7 = 0.1) and the coefficient ai =0.1 are such that the maximum 
rate of reaction is max{fi) = 0.025. In the case of ignition reaction the threshold was set up Cs — 0.6 while the 
exponent on the general Fisher non linearity is m = 10. Having fixed the grid step dx = 0.025 and the time step 
dt = 0.0005 we have run different simulations varying the coefficient Uq. 



By increasing Uq the front interface becomes more and more uneven. When Uq tends to one, it becomes 
discontinuous. Moreover, the baryccntric velocity of burnt fraction decreases by increasing Uo . 
These simulations were then repeated by solving equations (8) and (10) in the case of FKPP reaction 
(Fig. 3). We set up in both cases D = 7 = 0.1 and the initial condition c(.T,to) = c+(x,to) — (^{xito)- 
The numerical method can be said to be validated since the three models arc completely equivalent when 
Uq = 0, where the three equations degenerate in the case of incompressible flow. Moreover, increasing Uq 
but keeping it fairly low, there are no particular differences between the results obtained from equations 
(8) and (10). In this range of Uq, the hypothesis pD constant is not far from the hypothesis D constant. 
Nevertheless results are different from those obtained through the model represented by equation (10). 
Peaks of the concentration appear in proximity of the density peaks of the fluid. Increasing Uq, wc increase 
the compressibility and the concentration appears more and more peaked. Here, 6 has values higher than 
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the unity and hence confirming that in such areas the reaction rate is negative. For Uq tending to unity, the 
differences between the three equations are particularly marked. The diffusion model becomes particularly 
important since it is found, as expected, a relevant difference between the results from equations (8) (9). 
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Figures. Comparison between results from equations (10), (9) and (8) after a transient of t = 15 by increasing 
Uo- The reaction rate is an FKPP type, the initial condition c{x,to) = c^{x,to) = 9{x,to) with xq = 20.5 and 
ctq = 10, the grid step dx — 0.025 and the time step dt = 0.0005. For equations (10) and (9) we set up D = 7 = 0.1 
while for equation (8) pD = 7 = 0.1. 



We consider now the ignition reaction and we study under which conditon the reaction blows off. That 
is when the reaction stops before beeing completed. The phenomena can be studied by following in time 
the mass fraction of burnt gases which is defined, in normalized form, as: 

j^^ c(x,t)p{x)dx 

mb[t) = (15) 

Jo " P{x)dx 

If the reaction does not blow off, mb saturates to one. Otherwise, it stabilize to a lower value, which 
means that the fractions are no more reacting. 




(a) * (b) ' 

Figure4. On figure (a) the mass of burnt gases (irib) as function of time by increasing mb{t = 0) when Uo ~ 0.5. 
On figure (b) the mass of burnt gases {rrib) as function of time by increasing Uo and fixed mb{t = 0). The ignition 
rate of reaction with threshold Cs = 0.7 is considered, the parameter 02 = ^-^^'^ y2 = 1.111 and the diffusion 
coefficient pD — 0.4. 



We set up an initial condition as in (14) and we have varied ao, having fixed Uq — 0.5, until the blow-off 
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phenomena ocurr. Results are shown on Fig.4 which confirms that, in the case of relatively high diffusion 

and slow reactions, if the initial burnt fraction is too low the reaction can quench. Then, if wc fix the 
initial burnt fraction, Fig.4 shows that by increasing the compressibility of the flow field (increasing Uq) 
the blow-off phenomena can disappear. In paricular, there is a maxium value of Uo S [0.75, 0.8], for which 

the reaction qucnchs with more difficulty. 

We have then studied the possibility of blow-off when the m*^ order Fisher's reaction is applied. The 
approach has been the same, so that we firstly fixed Uq = 0.5 studying a range of ctq for which the blow 
off appears and than we fixed ctq studying the effects of compressibility. Results (see Fig. 5) shows that 
the rection never blows off. 




Figures. On subfigure (a) the mass of burnt gases {nib) as function of time by increasing mb{t = 0) when 
Uo = 0.5. On subfigure (b) the mass of burnt gases (mb) as function of time by increasing Uo and fixed mb{t = 0). 
The general mth order with m = 10 reaction is considered. The parameter as so that max[f2{c)] = max[f3{c)] 
and pD = 0.4 



4 Conclusion 

This brief analysis has helped us to shed light on some features of the models that are typically used 
in the study of population dynamics and in combustion in the form of an advection-reaction-diffusion 
equation. In particular, we showed that equation 6 applied to compressible flow fields leads to peaks of 
the concentration greater then unity. This means that the rate of reaction can take negative values. 
Since combustion is an irreversible process, a suitable reaction term has to be positive or zero and thus 
this model does not seem appropriate. We also studied the blow off phenomenon which is typical for the 
ignition like reactions. We have highlighted then the role of compressibility of the flow fleld. For slow and 
moderate compressible flows {Uq < 0.8), the efficiency of the reaction increases. Finally, we have studied 
the blow off phenomena when general mth order Fisher reaction is applied. For these kind of reactions 
this phenomenon does not occur. 
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